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Abstract 

The understanding of dynamics and functioning of biological membranes and in particular of 
membrane embedded proteins is one of the most fundamental problems and challenges in modern 
biology and biophysics. In particular the impact of membrane composition and properties and of 
structure and dynamics of the surrounding hydration water on protein function is an upcoming 
hot topic, which can be addressed by modern experimental and computational techniques. Corre- 
lated molecular motions might play a crucial role for the understanding of, for instance, transport 
processes and elastic properties, and might be relevant for protein function. Experimentally that 
involves determining dispersion relations for the different molecular components, i.e., the length 
scale dependent excitation frequencies and relaxation rates. Only very few experimental techniques 
can access dynamical properties in biological materials on the nanometer scale, and resolve dynam- 
ics of lipid molecules, hydration water molecules and proteins and the interaction between them. In 
this context, inelastic neutron scattering turned out to be a very powerful tool to study dynamics 
and interactions in biomolecular materials up to relevant nanosecond time scales and down to the 
nanometer length scale. We review and discuss inelastic neutron scattering experiments to study 
membrane elasticity and protein-protein interactions of membrane embedded proteins. 
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I. INTRODUCTION 



Dynamics in biological membranes occur on a wide range of length and time scales and 
involves interactions between the different constituents, such as lipids, cholesterol, peptides 
and proteins, as sketched in Figure [H Often, the structure of these systems is relatively 
well known and the corresponding experimental techniques, such as x-ray crystallography. 
Nuclear Magnetic Resonance (NMR), and also Atomic Force Microscopy (AFM) have devel- 
oped standard techniques in many disciplines. The access to dynamical properties turned 
out to be more difficult. Dynamical properties are often less well understood in biomolec- 
ular systems, but are important for many fundamental biomaterial properties such as, e.g., 
elasticity properties and interaction forces, and might determine or strongly affect certain 
functional aspects, such as diffusion and parallel and perpendicular transport through a bi- 
layer, and also protein function. New developments and improvements in neutron scattering 
instrumentation, sample preparation and environments and, eventually, the more and more 
powerful neutron sources open up the possibility to study collective molecular motions on 
lateral length scales between micrometers down to a few Angstroems (A). The fluctuations 
are quantified by measuring the corresponding dispersion relations, i.e., the wave vector- 
dependence of the excitation frequencies or relaxation rates. Because biological materials 
lack an overall crystal structure, in order to fully characterize the fluctuations and to com- 
pare experimental results with membrane theories, the measurements must cover a very 
large range of length and time scales The combination of various inelastic scattering 
techniques, such as inelastic neutron and x-ray scattering, and dynamic light scattering, en- 
larges the window of accessible momentum and energy transfers - or better: accessible length 
and time scales - and allows one to study structure and dynamics from the nearest-neighbor 
distances of lipid molecules to more than 100 nanometers (three orders of magnitude), cov- 
ering time scales from about one tenth of a picoseconds to almost 1 second (twelve orders 
of magnitude). The spectrum of fluctuations in model membranes for example covers the 
long wavelength undulation and bending modes of the bilayer with typical relaxation times 
of nanoseconds and lateral length scales of several hundred lipid molecules, down to the 
short-wavelength, picosecond density fluctuations involving neighboring lipid molecules. 

Even that membranes are now studied extensively for more than three decades the 
field is currently boosted by the inter disciplinary interest in life sciences. New and more 
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powerful techniques might give a quantitative access to molecular properties of the bilay- 
ers. Often, phospholipid bilayers are used as model systems to study fundamental properties 
3|, y, |5|. But the trend definitely goes to more complex and relevant systems, decorated with 
cholesterol, peptides and proteins to eventually develop a better understanding of biological 
membranes. Different model systems are used in experiments and Molecular Dynamics Sim- 
ulations (MD). In scattering experiments, the weak signals arising from dynamical modes 
can be multiplexed by using stacked membrane systems and large system sizes, i.e., mem- 
brane patches in the order of micrometers and larger. The advantage is the large sampling 
rate, which makes very small signals visible, and the large ensemble average, which averages 
over defects. The 'muli-scale' character of biological membranes, i.e., that relevant dynam- 
ics occur on a large range of length and time scales , poses particular problems 
because different techniques and probes must be applied to access the different lengths and 
frequencies. Scattering techniques are particularly powerful in spatially ordered systems and 
have been very successfully applied to study structure and dynamics (phonons) in crystals. 
In soft-matter and biology, a periodic structure is often missing or not well developed and 
new concepts and approaches are required to study structural and dynamical properties. 
In contrast to computer simulations, where different molecular components and dynamics 
can well be extracted and distinguished, in experiments there usually is a superposition 
of signals of different molecular components and interactions, and also single-molecule and 
collective dynamics. By using different sample geometries, these signals might be separated 
by their position in reciprocal and frequency space. Selective deuteration is used to dis- 
tinguish between incoherent and coherent dynamics. Because of more and more powerful 
spectrometers, smaller and smaller system dimensions and less and less material can be used 
in experiments. But experiments are still orders of magnitude away from single molecule 
spectroscopy of, e.g., membrane embedded proteins, which might be envisioned in the next 
decade. 

Even that computer power and access to powerful clusters has drastically increased in 
the last couple of years, MD simulations are still limited when it comes to large system sizes. 
Standard membrane systems are single bilayers, and contain several hundred lipid molecules, 
leading to patch sizes of about 100 A, and equilibrated for about 100 nanoseconds 
The system size, D, produces a cut off wavelength (A = 2t\: / D) for the dynamical modes 
that can be observed and quantified. To properly determine an excitation frequency or 
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relaxation rate from the time evolution of the scattering functions, the simulations must 
cover at least two orders of magnitude more in maximal accessible time, thereby limiting 
the accessible frequencies to about 1 ns. So simulations are still very much limited to fast 
short wavelength dynamics. When (single) proteins are simulated, the small system size 
prevents to study the influence of bilayer dynamics to protein dynamics and function, and 
also a possible protein-protein interaction, which might be relevant for biological function. 



Ana^ 
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^ical theory for membranes is usually based on the theory developed for liquid crystals 
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13| . It is challenging how microscopic properties, such as the lipid composition of 
a bilayer, affects macroscopic behavior. However, theories are very important to determine 
properties from fitting experimental and computational data. 

Collective molecular motions impact, e.g., on properties and functionalities of artificial 
and biological membranes, such as elasticity, transport processes, and inter-protein interac- 
tion. Two examples will be discussed in detail in the following: 

1. Mesoscopic shape fluctuations in aligned multi lamellar stacks of DMPC (1,2- 
dimyristoyl-sn-glycero-3-phoshatidylcholine) bilayers were studied using the neutron 
spin-echo technique. From the dispersion relation in the fluid phase, values for the bi- 
layer bending rigidity /t, the compressional modulus of the stacks B, and the effective 
sliding viscosity 773 could be determined. This technique offers a novel approach to 
quantify the elasticity parameters in membranes by direct measurement of dynamical 
properties and also the impact of collective molecular motions on membrane properties 



2. Very recently, interprotein motions in a carboxymyoglobin protein crystal were re- 
ported from a molecular dynamics simulation [Phys. Rev. Lett. 100, 138102 (2008)]. 
Experimental evidence for a cooperative long range protein-protein interaction in pur- 
ple membrane (PM) was found by inelastic neutron scattering. The dynamics was 
quantified by measuring the spectrum of the acoustic phonons in the 2d Bacteri- 
orhodopsin (BR) protein lattice. The data were compared to an analytical model and 
the effective spring constant for the interaction between protein trimers was deter- 
mined to be /c = 53.49 N/m. The experimental results are in very good agreement to 



the computer simulations, which reported an interaction energy of 1 meV 



Iq . In this 



case, inelastic neutron scattering was used to study interactions between constituents 
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of a biological membrane. 



Transport t 
layers [l6|, Il7 . 



iroiigh membranes is one of the most fundamental functionalities of the bi- 
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191]. While transport channels or transport proteins are responsible for 



the transport of larger molecules, small molecules, such as water, can pass the bilayer with 
a certain probability. Despite the large body of experimental and theoretical work, it turns 
out to be very difficult to model the permeation of small molecules through bilayers on 
a microscopic level because of the large number of parameters involved. The solubility- 
diffusion model and the transient pore model are possible mechanisms for permeation of 
small molecules [isl: Pores seem to be the dominant permeation mechanism for ions for 
thin membranes. Ion permeation by partitioning and diffusion seems to become of greater 
importance as membrane thickness increases. Neutral molecules, such as water, cross by the 
solubility-diffusion mechanism because of their high solubility in the hydrocarbon phase. 



Recently Nagle et al. [20( and Mathai et al. [21 1 presented in a combined study the the- 
oretical and experimental framework of a a three layer theory for the passive permeability 
through bilayers. They link permeability to structural properties and strongly correlate it 
with the area per lipid^, rather than with other structural quantities such as the thickness. 
Decades ago Trauble [2^ (and later addressed by Marsh [23|) discussed the movement of 
molecules across membranes in terms of thermal fluctuations in the hydrocarbon chains of 
the membrane lipids. The thermal motion of the hydrocarbon chains results in the formation 
of conformational isomers, so-called kink-isomers of the hydrocarbon chains. "Kinks" may 
be pictured as mobile structural defects which represent small, mobile free volumes in the 
hydrocarbon phase of the membrane. The free volume needed for kink creation and move- 
ment might be created by propagating in-plane density fluctuations of the lipid acyl chains 



24| . In a first inelastic neutron scattering experiment, the corresponding dispersion relation 



of the acyl tail dynamics in DMPC bilayers showed a drastic softening in the presence of 
ethanol molecules, which are known to increase bilayer permeability. Future experiments 
will elucidate the infiuence of ethanol on the collective dynamics in phospholipid bilayers. 

II. COLLECTIVE VS. SINGLE-MOLECULE DYNAMICS 



The spectrum of fiuctuations in biomimetic and bi olog ical membranes covers a large range 



of time and length scales 
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28|, 



29|, 



30|, ranging from the long wavelength 



undulation and bending modes of the bilayer with typical relaxation times of nanoseconds 
and lateral length scales of several hundreds lipid molecules to the short wavelength density 
fluctuations in the picosecond range on nearest neighbor distances of lipid molecules. Local 
dynamics in lipid bilayers, i.e., dynamics of individual lipid molecules as vibration, rota- 
tion, libration (hindered rotation) and diffusion, have been investigated by, e.g., incoherent 



neutron scattering 
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27, 
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291] and nuclear magnetic resonance 3l|, 12^ to determine 



the short wavelength translational and rotational diffusion constant. Collective undulation 



modes have been investigatec 



dynamical light scattering |34l. l35l. l36|. 



using 
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3311 and 



Atomic and molecular motions in membranes can be classified as local, autocorrelated, 
and collective, pair correlated dynamics. Within the scattering formalism, autocorrelated 
dynamics is described by the incoherent scattering function, SincoherentiQ,^^), while the co- 
herent scattering, Scoherentil,^), describes dynamics involving different molecules or atoms. 
Figure [2] exemplary depicts some of the local and collective modes in a phospholipid bilayer. 
Rotational and lateral diffusion, vibrations and rotations of the single lipid molecules can 
be investigated by e.g., incoherent inelastic neutron scattering, nuclear magnetic resonance 
or dielectric spectroscopy. On the other hand only coherent inelastic neutron scattering or 
inelastic x-ray scattering are able to elucidate the collective excitations such as, e.g., the 
short wavelength density fluctuations or undulation modes of the bilayer. 

Experimentally, selective deuteration is used to emphasize the incoherent, respective co- 
herent scattering over other contributions to the total sample scattering. While in proto- 
nated samples the incoherent scattering is normally dominant and the time-autocorrelation 
function of individual scatterers is accessible in neutron scattering experiments, (partial) 
deuteration emphasizes the coherent scattering and gives access to collective motions by 
probing the pair correlation function. Computer simulations offer direct access to incoher- 
ent and coherent properties by calculation of SincoherentiQ, ^) and ScoherentiQ, t^) from the time 
evolution of atomic and molecular coordinates. Figure [3] visualizes some of the length and 
time scales involved for collective dynamics and the corresponding functional aspects. While 
permeability of bilayers occurs on distances of neighboring lipid molecules, elasticity might 
involve hundreds of membrane molecules. Coherent intra-protein dynamics will most likely 
be faster than inter-protein dynamics, which involves larger distances and in most cases a 
lipid mediated, elastic interaction. 
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III. EXPERIMENTAL 



Only recently, the first inelastic scattering experiments in phospholipid bilayers to de- 
termine collective motions of the lipid acyl chains and in particular the short wavelength 
dispersion relation have been performed using inelastic X-ray js^] and neutron 24] scatter- 
ing techniques. Note that only scattering experiments give wave vector resolved access to 
dynamical properties, what is important to associate excitation frequencies and relaxation 
times with specific molecular components and motions. In the case of single membranes the 
scattering signal is usually not sufficient for a quantitative study of the inelastic scattering. 
To maximize the scattering signal, multi lamellar samples composed of stacks of several 
thousands of lipid bilayers separated by layers of water, resulting in a structure of smectic 
A symmetry, were prepared. The high orientational order of the samples which gives rise to 
pronounced Bragg peaks and excitations is a prerequisite to a proper analysis of the corre- 
sponding correlation functions. Highly oriented multi lamellar membrane stacks of several 
thousands of lipid bilayers were prepared by spreading lipid solution of typically 25 mg/ml 
lipid in trifluoroethylene/chloroform (1:1) on 2" silicon wafers. About twenty such wafers 
separated by small air gaps were combined and aligned with respect to each other to create 
a "sandwich sample" consisting of several thousands of highly oriented lipid bilayers with a 
total mosaicity about 0.5°), and a total mass of about 400 mg of deuterated DMPC The 
mosaicity of the sandwich is composed of the alignment of the bilayers within the stack on 
one wafer and of the orientation of different wafers with respect to each other. The use of well 
oriented samples leads to well localized elastic and inelastic signal in reciprocal space and 
allows to distinguish motions in the plane of the membranes (g||) and perpendicular to the 
bilayers (g^). During the experiments, the membranes were kept in a "Humidity Chamber" 
to control temperature and humidity and hydrated with D2O from the vapor phase. The 
collective motions of the lipid acyl chains have been emphasizes over other contributions to 
the inelastic scattering signal by using partially, chain deuterated lipids (DMPC -d54). The 
experiments that will be discussed here, were conducted on two different types of neutron 
spectrometers, namely triple-axis and spin-echo spectrometers. The accessible length and 
time scales of these spectrometers are discussed in more detail in Ref. The present paper 
focuses on the dynamics-property and dynamics-function aspects of the experiments. 
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A. Spin-echo Spectrometry 



Fluctuations on the mesoscopic scale are determined by the elasticity parameters of the 
bilayers, i.e., the compressibility of the stacked membranes, B, and the bending modulus 
n. The relaxations in this regime are in the nano-second time-range with accompanying 
small g- vales. Spin-echo spectrometers turned out to be highly suited for these experiments. 
The spin-echo technique offers extremely high energy resolution from Larmor tagging the 
neutrons. A neutron spin echo measurement is in essence a measurement of neutron polar- 



ization |38|. A polarized neutron beam passes through a magnetic field perpendicular to the 
neutron polarization. The neutron spin precesses before arriving at the sample, acquiring a 
precession angle ipi. At the sample, the beam is scattered before passing through a second 
arm, acquiring an additional precession angle in the reversed sense. For elastic scattering 
the total precession angle is Aip = y^i — </?2 = for all incoming neutron velocities. If the 
neutron scatters inelastically by a small energy transfer hu, there will be a linear change 
= T ■ uj with r being a real time in the case of quasielastic scattering. The spin-echo 
technique thus works in the time domain and measures the intermediate scattering function 
S(g||,t) in contrast to the three-axis technique. For a quasielastic response, assumed to have 
Lorentzian lineshape with half- width F, the polarization will then show a single exponential 
decay PNSE = P.e^*. 



B. Triple-axis spectrometry 

Fast motions in the ps time range due to sound propagation in the plane of the bilayer 
are best measured on triple-axis spectrometers. The energy of the incident and scattered 
neutrons is determined by Bragg scattering from crystal monochromators, (graphite in most 
cases). Advantages of triple-axis spectrometers are their relatively simple design and op- 
eration and the efficient use of the incoming neutron flux to the examination of particular 
points in (g, cu) space. By varying the three axes of the instrument, the axes of rotation of 
the monochromator, the sample and the analyzer, the wave vectors ki and kf and the ener- 
gies Ej and Ej of the incident and the scattered neutrons, respectively, can be determined. 
The momentum transfer to the sample, and the energy transfer, Tiuj, are then defined by the 
laws of momentum and energy conservation to q = kf — ki and /ia;=Ej-E/. The accessible 



8 



(g, uj) range is just limited by the range of incident neutron energies offered by the neutron 
guide as well as by mechanical restrictions of the spectrometer. 



IV. ELASTIC PROPERTIES 

According to linear smectic elasticity theory jsQi. I40l| thermal fluctuations in the fluid phase 
of the membrane are governed by the free energy functional (Hamiltonian) 39|, |40|, |4l| : 

H = j/r X: - u^f + -K (vSn„) j , (l) 

where k denotes the bilayer bending rigidity, A the area in the xy-plane, the number 
of bilayers, and the deviation from the average position n d oi the n-th bilayer, d is 
the lamellar spacing. B and K = njd are elastic coefficients, governing the compressional 
and bending modes of the smectic phase, respectively. A fundamental length scale in these 



systems is given by the smectic penetration length A = ■JKJB. Aligned lipid bilayers allow 
a separate determination of both parameters K and B 42l. |43 1 . 

The spin-echo experiments were carried out at the INll and IN15 spectrometers, at 
the cold source of the high flux reactor of the Institut Laue-Langevin (ILL) in Grenoble, 
France. Wavelength bands centered at A=7.4A and A=14A with AA/A ~ 0.15 (FWHM), 
respectively, have been set by a velocity selector. The intermediate scattering function 
S(q||,t) was measured for spin-echo times of 0.001ns<t<20ns for INll and 0.01ns<t<200ns 
for IN15. Data have been taken at three different temperatures, at 19°C, in the gel (ripple, 
P/3') phase of the phospholipid bilayers, at 22°C, just above the temperature of the main 
transition in deuterated DMPC-d54 (at T^ ~ 21.5°C), and at 30°C, far in the fluid L^ phase 
of the membranes and above the regime of so-called anomalous swelling. The corresponding 
lamellar d spacings were (i=56 A, 60 A and 54 A (gel, 22°C and fluid), respectively. Two 
relaxation processes, one at about 10 ns (ri) and a second, slower process at about 100 ns 
{T2) were observed. The relaxation rates rf ^ and in the gel and the fluid phase are 
depicted in Fig. H] (a) and (b). Both relaxation branches are dispersive. The fast process 
shows a increase at small q|| values and a bend at about q|| R:i0.015A~^. The dispersion in 
the gel phase and close to the phase transition in Fig. H] (b) appear to be more pronounced as 
compared to 30°C dispersion. A soft mode appeared in the T=22°C dispersion, indicating 
a signiflcant softening of the bilayer at a well defined wave number. The slow branches at 
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45|, 



T=19°C and 30°C also show increasing relaxation rates with increasing q|| values, but with 
a distinct non-polynomial behavior. 

The dispersion relation of the fast branch with relaxation rates between 1 and 10 ns can be 
attributed to undulation dynamics. Qualitatively, at very small g||-values, the membranes 
behave as liquid films and their dynamics is basically determined by the viscosity of the 
water layer in between the stacked membranes {film regime). With increasing q\\ there is a 
transition into a bulk- elasticity regime where the dynamics depends on the elastic properties 
of the lipid bilayers. At this point, the dispersion bifurcates in two relaxation branches. 
The faster one, which is out the experimentally accessible time window of the spin-echo 
spectrometer, is mainly determined by the compressional modules B. The slower one, which 
is observed here, can be assigned to the bending modulus k. The slow dispersion branch wit 
relaxation rates of about 100 ns could be attributed to a surface relaxation mode ; 
which is particular to a stack of membranes. 

Following the idea of Ribotta ISj, the relaxation rates of the undulations can be described 

(?73 is the layer sliding viscosity), and following results were obtained: k, = 14.8 ± Sk^T, 
A = 10.3 ±2.3 A, r/3 = 0.016 ±0.0006 Pa s. B was calculated to B=1.08 lO^m^ (rf=54 A). 
These values agree quite well with values reported in the literature 

Q, y, The 

resulting effective sliding viscosity of the membrane system 773 was found to be 16 times 
higher than that of water, what points to different properties of the interstitial hydration 
water, as compared to bulk water. Note that Eq. [2] does not describe a pure undulation 
mode, which is probed at q^ values of = 27r/d, only. If the scattering is probed at finite 
components 6qz '■= {qz — ^n/d) or measured with a relaxed q^ resolution, there is a mixing 
of baroclinic modes, which are distinctly slower than a pure undulation because they involve 
a relocation motion of the water layer. The parameter D describes an effective finite-size 
cutoff- length, which was related to the instrumental resolution D = vr/Ag^. 

From the soft-mode in the fluid dispersion, it seems that the well known softening of 
phospholipid membrane upon approaching the main phase transition temperature from the 
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occurs on 



fluid phase, i.e., the regime of "critical swelling" or "anomalous swelling" 
a well defined length scale, only (27r/g|| ^ 420 A). Using atomic force microsopy (AFM), 
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pronounced ripples were observed in the P^/ phase of stacked membranes, with an increasing 
ripple per iodicity when approaching the temperature of the main transition at Tm = 
24 °C 49 1 . Close to T^, coexisting metastable ripples with 2Ar ^ 420 A were observed. 



The experiments therefore point to coexisting nanodomains with sizes of less than 50 nm 
in the range of critical swelling 
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52l | of phospholipid bilayers. The existence of 



coexisting small gel and fluid domains has also been argued by preceding AFM investigations 
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54j to compensate the large stress which occurs at Tm due to the volume difference of 



the two phases. The softening in the range of the phase transition is most likely coupled 
to the occurrence of these metastable 2A,. ripples. Bending of the bilayers might occur 
mainly in the interfaces between two metastable ripples where the bending modulus can 
be expected to be softer because structure and interactions are likely to be much less well 
defined in the interface. Not much energy would be needed to slightly change the tilt angle 



between to ripple flanks. Note that the relation between critical swelling and so: 
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tening of the 



55| and the 



bilayers, and the formation of a low temperature ripple phase is not trivial 
anomalous swelling is not directly coupled to the formation of ripples. Inelastic experiments 
in different systems are definitely needed to better understand the impact of the soft-mode 
in the DMPC dispersion curve. The most important contribution of the NSE technique 
is that the 'stiffness' of the bilayers is determined as a function of internal length scale. It 
seems that there is no global softening, but softening occurs on of certain length scales, only. 

V. COLLECTIVE PROTEIN INTERACTION 



The high protein concentration in biological membranes might lead to long-range protein- 
protein interactions, on which there have been speculations, already some time ago j^. 
Recently, interprotein motions in a car 



molecular dynamics simulation 
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:boxymyoglobin protein crystal were reported from a 
571]. Motions in proteins occur on various length and 
time scales p, |7|], and the functional behavior of membrane proteins is likely to depend on 
the lipid bilayer composition and physical properties, such as hydrophobic thickness and 
elastic moduli. How the variety of inter- and intra-protein motions, occuring over different 
time and length scales, interact to result in a functioning biological system remains an open 
field for those working at the interface of physics and biology. Purple Membrane (PM) 
occurs naturally in the form of a two-dimensional crystal, consisting of 75% (wt/wt) of a 
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single protein, Bacteriorhodopsin (BR), that functions as a light-activated proton pump, 
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BR is a proton 



and 25% various lipid species (mostly phospho- and glyco-lipids) 
transporting membrane protein, formed of seven trans-membrane alpha-helices arranged 
around the photosensitive retinal molecule. The protein in the lipid matrix is organized in 
trimers that form a highly ordered 2d hexagonal lattice with lattice parameter a ~ 62 A, as 
depicted in Fig. [5] (a). 

The experiments were performed on the IN12 cold-triple-axis spectrometer at the Institut 
Laue Langevin (Grenoble, France). It allows the measurement of diffraction and inelastic 
scattering in the same run without changing the set-up, which is crucial to assign dynamical 
modes to structural properties and molecular components. Correlations and motions in 
membranes are often well separated in reciprocal space because of the largely different length 
and time scales involved. The prominent distances in PM, such as lipid-lipid and BR- 
BR monomers and trimers for instance lead to spatially well separated signals. The same 
holds for the different time scales involved from the picosecond (molecular reorientations) 
to the nano- or microsecond (membrane undulations, large protein motions). The use of 
oriented samples further allows to separate correlations in the plane of the membranes, and 
perpendicular to the bilayers. Dynamics between different protein trimers is expected to 
be dominant where the 2d BR diffraction pattern is observed, i.e., in a g|| range of about 
0.1 to 0.6 A-\ 

The excitation spectrum of the 2d protein lattice was modeled analytically, taking the 
protein trimers as the centers of a primitive hexagonal lattice with lattice constant a = 
62 A. The model is depicted in Fig. \5la). The basic hexagonal translations are marked by 
arrows. The interaction between the protein trimers is contained in springs with an effective 
(longitudinal) spring constant k (Fig. [5t^b)). The calculated Ci{q,uj) is shown in Fig. [6](b). 
The statistical average leads to a superposition of the different phonon branches, which start 
and end in the hexagonal Bragg peaks (at huj = 0). The absolute phonon energies can not 
be determined from the model, but depend on the coupling constant k. So the energy of 
the phonon curves in Fig. E^b) was scaled to match the experiment. Note that because 
the proteins trimers were treated as dots with an effective mass of Mtr, the calculation 
does not include any contributions from intra-protein or intra-trimer dynamics, i.e., possible 
optical modes and phonons. The longitudinal phonon spectrum Ci{q,uj) = {u'^ / q^)S{q,uj) 
for q\\ values between 0.34 A~^ and 0.46 A~^ is shown in FigureE](a). The most pronounced 
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phonon branches from Figure [6](b) were plotted as solid lines in part (a) for comparison. 
In the q\\ range between 0.34 and about 0.43 A"^, the experiment well reproduces the 
calculated phonon curves. The agreement is less good for q\\ values above 0.43 A~^. Only 
the strongest phonon branches are visible in the data, the weaker branches can most likely 
not be resolved from the background. 

Experiment and calculation can be compared to the MD simulations in Ref. 561], where 
a peak in the energy spectrum at 1 meV was identified as a translational intermolecular 
protein:protein interaction vibration in a carboxymyoglobin protein crystal. The energy 
agrees well with the energy of the zone boundary phonon in PM of 1.02 meV, as shown in 
Fig. O The good agreement of the energy values in the two systems most likely stems from 
the very high protein density in PM, which makes it almost crystal like. The computational 
work thus strongly supports the interpretation of our data as collective protein:protein ex- 
citations. The commonly assumed interaction mechanism between inclusions in membranes 



is a lipid-mediated interaction due to local distortions of the lipid bilayer 59|, l60|, l6l|, l62|, |63| , 



with a strong dependence on the bilayer properties, in particular elastic properties. The 
PM might, however, be a special case because there are very few lipids between neighboring 



BR proteins 



6J]. While the nature of the interaction still will be mainly elastic, it is not 



likely to be purely lipid-mediated but for the most part a direct protein-protein interaction. 
The strength of the interaction can be determined from the data in Fig. [61 The energy 
of the zone-boundary phonon at the M-point of the hexagonal Brillouin zone (for instance 
at a g|| value of 0.35 A~^) relates to the coupling constant by MtrU}"^ = Qk. Because this 
energy is determined as hu = 1.02 meV, the effective protein-protein spring constant k is 
calculated to k = 53.49 N/m 66||. On the microscopic level, displacing the BR trimer by 1 
A yields a force between neighboring trimers of 5.3 nN. There is therefore strong protein- 
protein communication in PM. Using the same approach, the spring constant for graphite 
for comparison is calculated to 27,000 N/m for the in-plane interaction, and 3.5 N/m for 
out-of-plane interactions. The force constant that we measure in PM thus is 1-2 orders of 
magnitude larger than the effective van-der-Waals force constant in graphite, but 2-3 orders 
of magnitude weaker than a C-C bond. 
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VI. CONCLUSION 



Inelastic scattering experiments give access to molecular dynamics and correlations in 
membranes. The experiments prove that relevant properties and functional aspects of the 
bilayers can be quantified from coherent dynamics. Future experiments will allow to study 
intra and inter protein dynamics of membrane embedded proteins and how bilayer compo- 
sition and physical properties affect protein function. Membranes with specific properties 
might be tailored for specific applications and protein function of membrane embedded 
proteins might be enhanced (or suppressed) by adapting the properties of the surrounding 
bilayer. 
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Figure Legends 

Figure [I] 

(Color online). Cartoon of a cell membrane. Courtesy of Mariana Ruiz Villarreal 
( http: / / commons.wikimedia.org/wiki /Image:CelLmembrane_detailed_diagram.svg ) . 



Figure 

(Color online). Some elementary dynamical modes in lipid bilayers. (a) Local modes 
include diffusion and vibrations, rotations and librations (hindered rotations) of single lipid 
molecules, (b) Collective excitations are coherent motions of several membrane molecules, 
such as the long wavelength undulation and bending modes of the membranes. 

Figure 

(Color online). Length and time scale for some collective motions and corresponding func- 
tional aspects. While permeability occurs on nearest neighbor distances of lipid molecules, 
membrane elasticity involves several hundreds of lipid molecules. Besides collective intra- 
protein dynamics, there might also be an interaction between different proteins in biological 
membranes, i.e., an inter-protein dynamics. 



Figure 



(Color online), (a) Dispersion relations at T=30°C. The solid line is a fit to Eq. ([2]). (b) 
Dispersion relations in the gel (19°C) and in the fluid phase (22°C). A pronounced soft mode 
is observed at qn_~0.015A~^ at 22°C (dotted vertical line). Solid lines in (b) are guides to 
the eye. (From 14|.) 



Figure 

(Color online), (a) BR trimers are arranged on a hexagonal lattice of lattice constant 
a ~ 62 A. (b) The interaction between the protein trimers is depicted as springs with 
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effective spring constant k. (From [1 



Figure 

(Color online), (a) Experimental Ci{q,uj), measured in the the third Brillouin zone, (b) 
Calculated excitation spectrum Ci{q,uj) in the range of the experimental data. The most 
strongest phonon branches are also shown in part (a) as solid lines. (From 15|.) 
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